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Abstract 

"Weak measurements" - involving a weak unitary interaction between a quantum system and 
a meter followed by a projective measurement - are investigated when the system has a non- 
Hermitian Hamiltonian. We show in particular how the standard definition of the "weak value" of 
an observable must be modified. These studies are undertaken in the context of bound state scat- 
tering theory, a non-Hermitian formalism for which the Hilbert spaces involved are unambiguously 
defined and the metric operators can be explicitly computed. Numerical examples are given for a 
model system. 
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I. INTRODUCTION 



The standard formulation of quantum mechanics requires physical observables to be 
mathematically given in terms of Hermitian operators. In the last decade theories with 
a non- Hermitian Hamiltonian have been extensively investigated jl]. The initial momen- 
tum was given by work concerning PT-symmetric Schrodinger operators It was hoped 
that the PT-symmetric Hamiltonians, which are complex but nevertheless possess a real 
spectrum, would provide an extension of standard quantum mechanics. It was later ar- 
gued however that these non-Hermitian operators could be mapped to Hermitian ones by a 
similarity transform Nevertheless the non-Hermitian framework remains useful. Indeed, 
from a fundamental perspective it opens up the possibility of doing quantum mechanics 
with non-standard inner products. This has practical consequences because many physical 
systems are naturally formulated in non-Hermitian terms 4[]. 

Scattering systems involving bound states in real potentials are such a case. While scat- 
tering problems in complex potentials have been prominent in the studies of non-Hermitian 
Hamiltonians,bound state scattering in short-range real potentials have been scarcely in- 
vestigated [5H7|. In this case, non-hermiticity arises from the boundary conditions imposed 
on the scattering functions, ultimately linked to the fact that the scattering solutions are 
not the full eigenstates of the exact Hamiltonian. From the physical point of view this exact 
Hamiltonian exists, but its eigenstates are unknown in practice, while the bound scatter- 
ing solutions are eigenstates of an effective Hamiltonian that is not Hermitian relative to 
the standard inner product. Hence in principle one should employ a biorthogonal basis, or 
equivalently obtain the metric operator in order to define the inner product relative to which 
the effective Hamiltonian becomes self-adjoint. 

In this paper, we will focus on the interaction between a bound state scattering system 
and a measurement device in a scheme popularly known as "Weak measurements" (WM). 



WM, introduced more than 2 decades ago 



have been receiving an increased attention 



these last 5 years, in particular as theoretical but also experimental tool aimed at inves- 
tigating fundamental problems in quantum mechanics (see 9| and Refs. therein). WM 
actually involve two steps: the first step is a weak interaction between the system and a 
"weak meter", the overall evolution being unitary. The second step is a standard projective 
measurement in which the system (at that point entangled with the weak meter) interacts 



2 



with a different measurement device. The state of the system is projected to a final post- 
measurement state, while the weak meter has picked up a phase depending on a quantity 
known as the "weak value" of the weakly measured observable. In the conjugate variable 
of the pointer, the phase-shift appears as as a shift in the probability distribution. This 
shift can be experimentally measured by obtaining the probability distribution of the weak 
pointer. 

The main issue when considering a weak measurement of a non-Hermitian system lies in 
the treatment of the coupling between the system and the weak meter. Indeed in a standard 
projective measurement the observed quantity is an eigenvalue, which is a real quantity not 
depending on the definition of the inner product. In a weak measurement, the observed 
quantity is a shift in the pointer proportional to the weak value, which, as will be seen 
below, is a renormalized transition element. It is therefore crucial in order to determine the 
weak value, to properly define the inner product and the physical Hilbert space. 

We will first briefly introduce weak measurements and give the usual formula for comput- 
ing weak values (Sec 2). We will then derive the weak value for non-Hermitian systems. In 
order to provide an unambiguous physical basis, this derivation will be done in the context 
of bound state scattering theory. We will thus explain why systems described with this for- 
malism are non-Hermitian in the 'physical' Hilbert space, leading to the definition of a new 
inner product and its associated Hilbert space, in which the weak values must be defined 
(Sec 3). We will then give in Sec 4 examples of weak value computations for a model bound 
scattering system. We will see that non-Hermitian issues must be incorporated explicitly in 
order to account for the correct shift in the weak measurement apparatus. Our concluding 
remarks will be given in Sec 5. 

II. WEAK MEASUREMENTS 

A standard quantum measurement, often represented by the projection of a premea- 
surement state of the system to an eigenstate of the measured system observable, actually 
involves a two-step procedure. First a unitary interaction between the measured system and 
the measurement apparatus results in a system-apparatus state entangled in the pointer 
basis. Then in a second step the entangled state is projected to a final post-measurement 
state correlating a unique pointer state with an eigenstate of the measured observable. 
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A weak measurement of an observable A proceeds differently. First a weak unitary 
interaction takes place between the system and the "weak" meter. The weakness of the 
interaction results in an entanglement [ji]] in which the different pointer states are nearly 
identical. Then a standard quantum measurement of a different observable takes place, 
resulting in the usual projection to an eigenstate of this second observable. Since the system 
and the weak apparatus were still entangled, the projection to a final state of the system 
also determines the quantum state of the weak meter. 

Rather than solving for the weak interaction in terms of the entangled states in the pointer 
basis, the standard approach 8| to weak measurements starts from a first order expansion 
of the interaction Hamiltonian. Let \ip(ti)) and be the initial states of the system 

and weak meter respectively just before they interact, and let us assume an interaction 
Hamiltonian of the form I(t) = f{t)AX where the system observable A is coupled to the 
weak pointer's position variable along the x axis (ie X \X) = X \X) for the pointer). f(t) is a 
smooth function of t vanishing outside the interval t_ < t < t + during which the interaction 
takes place and obeying J t * + f(t)dt = g where g is the mean effective coupling strength. 
Neglecting the self-evolution of the system and meter during the time interval t + — t_, the 
unitary evolution generated by I(t) brings the initial state = IV'(^)) I^CO) ^° 

m + )) = e-^ AX MU))MU)). (1) 

A projective measurement of another observable B of the system is made immediately after. 
The system state is projected to one of the eigenstates of B; among the possible outcomes 
we select only the cases in which the final state is \/3f) . The standard approach consists in 
expanding the exponential to first order in g and then compute the projection 

(X\ <J) f \ ~ ^)> (X\ exp (-ig M^M xj Mtt )} (2) 

The term 

is known as the weak value of A. Eq. (j2J) indicates that the weak meter has picked up a 
phase (in configuration space), or alternatively a shift (in momentum space) proportional to 
the weak value of A, given the initial (known as "preselected") state, and the final ("post- 
selected" ) state obtained after having made a standard measurement of another observable 
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B. Note that the weak value can be a complex number, implying different shifts can be 



111- 



observed in the conjugate variables of the meter 

The derivation of Eq. (j2J) involves several approximations (see eg [10|) that will not be 
discussed here. A necessary (but not sufficient) condition is that g(A) w is small, gener- 
ally implying that the coupling g must be vanishingly small (because (A) w is generally 
large). The important point, from a physical perspective, is that measuring the weak meter 
wavefunction allows to obtain information, encoded in (A) w , on the system observable A 
without making a full quantum measurement of that observable. Instead another, possibly 
incompatible property B is measured. 

The applications and interpretations of WM are out of the scope of this work. Our focus 
here lies in the weak measurement of a system described in a non-hermitian framework. 
Indeed, accounting for WM involves treating - though only to first order - the interaction 
between the non-hermitian system and a measurement device. Contrary to a standard 
measurement, in which case the outcome would be an eigenvalue, the measurement device is 
shifted by the weak value. However the definition ([3]) of the weak value is valid in standard 
(Hermitian) quantum mechanics. For a system described in a quasi-Hermitian framework 
(A) w must be computed in the correct Hilbert space, endowed with a non-standard inner 
product, as will be seen below. 

III. NON-HERMITIAN FORMALISM 
A. General remarks 

The weak value as given by Eq ([3]) needs to be modified for systems described in a 
non-Hermitian setting. The rationale, well-known to practitioners of PT-symmetric/quasi- 
Hermitian quantum mechanics, is that the inner product needs to be replaced. Given the 



controversies surrounding the physical interpretation of non-Hermitian systems [12|, our 
approach will consist in working with a system - or rather family of systems, those involving 
bound state scattering - that has an important advantage: the non-Hermitian aspects appear 
because one is led to work with wavefunctions defined on a modified configuration space. 
This means that while the scattering Hamiltonian is non-Hermitian, there is in principle 
an underlying exact Hamiltonian (though untractable in practice). As a result there is no 
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ambiguity when delaing with the conceptual aspects surrounding non-Hermiticity. 

We will therefore first give a brief presentation of bound state scattering, exposed previ- 
ously in Ref js| , and discuss its non-Hermitian aspects in order to derive the formula for the 
weak value in non-Hermitian systems, given by Eq below. 



B. Bound-state scattering 

For definiteness, let us consider 2 particles, a light particle and a massive compound tar- 
get, attracted by a long-range radial field. The scattering between the particles is described 
by a short-range potential. Letting H e denote the exact Hamiltonian in the center of mass, 
we assume H e can be split as 

H e = H + V (4) 

where Hq is the Hamiltonian of the light particle in the long-rang field while V contains all 
the short range interactions between the light particle and the target. We further assume 
that 

(r'\V\r) = 9(r -r')V9(r -r), (5) 

ie V vanishes outside some small radius ro (9 is the step function). The total energy E can 
be partitioned as 

E = e l + e i (6) 

where £j is the internal energy of the target (depending on the target quantum state) and 
6i is the energy of the light particle. The eigenstates of Ho are given by 

^{E)) = |/*fe)> |i(e,)> ; (7) 

fi(ei,r) = (r |/i(ei)) is the eigenfunction of the radial part of Hq whereas the 'target' state 
|i(£j)) includes all the other degrees of freedom, including the non-radial ones of the colliding 
particle (a handy notation given that the angular momenta of the particles are usually 
coupled). The target states are orthogonal, (i\ j) = 5ij. For bound states fi{e^r) vanishes 
at and +oo (whenever E is an eigenvalue of H Q ). 

The label % defines the scattering channel. In each channel the standing-wave solutions 
are given by the Lippmann-Schwinger equations of scattering theory as 

mE)) = \HE)) + G (E)K(E) ME)) (8) 
6 



where Gq(E) is the principal- value Green's function and K the reaction (scattering) operator 
for standing waves linked to the familiar S matrix by a Cayley transform [jjj . The difference 



here with standard scattering theory is that the bound channels are included explicitly [20]. 
Both (r| <fri(E)) and (r\ ip?(E)) diverge as r — > oo for an arbitrary value of E. A bound state 
appears when the superposition 

mE)) = Y,z l {E)mE)) (9) 

i 

converges as r — > oo. This happens for discrete values of the energy obtained, along with 
the expansion coefficients Zi(E), by imposing the boundary conditions. 

While H e is undoubtedly Hermitian relative to the standard inner product, its eigenfunc- 
tions cannot be computed from Eqs ([H])-© because the formal expansion of Go over the 
eigenstates of Hq is intractable. Instead the scattering formulation consists in obtaining a 
closed form expression of Go but valid only outside the reaction zone, ie for r > r$. Indeed 
from the scattering viewpoint, whatever happens within the reaction zone is encoded in the 
phase-shifts. The wavefunction ([9]) outside the reaction zone takes the form 



|V(£)) = 5>(£) 



fi(e i} r) \i) + ^2g j (e jl r) \j) K j{ 



r > r (10) 

j 

where are the on-shell elements of the scattering matrix, which are assumed to be known. 
g(r) is, like f(r) introduced in eq (GO), a solution of the radial part of Hq but it is irregular 
at the origin. 

The scattering state \ip{E)) of (TTOT) is the part for r > r of the exact solution \ijj G {E)), 
and not an approximation to it. But within the scattering formulation the 'inner' part of 
\ip e (E)) for r < r does not exist: all meaningful quantities are defined radially on [r ,oo[. 
As a consequence, 

^(E 2 )) = 5 ElE2 + fi ElE2 (l - 5 ElE2 ) (11) 
ie the scattering states are normalized to 1, but are not orthogonal; this is due to the fact 



that the boundary conditions at r = r are not identical for all the \ip(E)) 
scattering states cannot be eigenstates of the Hermitian operator 



Hence the 



H = Y,E\1>m(4>W\ (12) 

E 

since H \ip(E)) ^ E \ip(E)) . A non-Hermitian Hamiltonian needs to be introduced instead. 



C. Non-Hermit ian aspects: metric, Hilbert spaces and operators 



Let 1-iph be the Hilbert space of standard quantum mechanics. Physical states are repre- 
sented by vectors in T-L p h- From a practical viewpoint, we may consider that the phase-shifts 
(or the K matrix elements) are known and the problem concerns the expansion of physical 
states in terms of the scattering solutions \ip{E)) . 

To this end let us introduce a non-Hermitian operator H and state vectors 
that 

H\i;(E)) = E\i;(E)) 



H- 



E 



m) we')) 



>EE> 



ip{E)j such 

(13) 
(14) 

(15) 



where { i/j(E)j , \ip(E))} forms a biorthogonal basis. It follows that we can write the follow- 
ing expansions: 

H = J2 E l^)> ($( E )\ H + = Y,E \^{E)) ^{E)\ . (16) 

E E 

H and H + are further linked by 

h = gH+g- 1 (17) 

where g is a Hermitian operator given by 

= £ mm q- 1 = £ |^)) (^)| • ( 18 ) 

E E 

Eq m is the defining relation of ouaei-Herooitieity Q Q being inwrtible and positive- 
definite |5J. 

This allows to define a Hilbert space "H endowed with a new inner product depending on 
the metric 

WE^iE^g = (HEi)\g^ \i>(E 2 )) = (i){Ei) \^{E 2 )) = 5 ElE2 . (19) 

Eq ( TTT|) indicates that H is Hermitian relative to this new inner product. Completeness of 
the biorthogonal basis allows to expand an arbitrary state of 7~l v h in terms of the \i/j(E)), ie 
the eigenstates of H span the entire Hilbert space of admissible physical states even if they 
do not form an orthogonal basis in 1-L p h 21] 
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Calculations involving the scattering states have to be performed in % rather than in % p u- 
Indeed although a physical state \oc k ) h is known in T-L p h, its expansion over the scattering 
eigenstates \a k ) = Yl,E a k(E) \ip{E)) is defined in % with the expansion coefficients given 
through 

a k (E) = 1^{E) \a k ) = (ip(E),a k ) g = (iP(E) \a k ) (20) 

where we have put \a k ) = Q~ l \a k ) . Note that in the underlying exact problem, there is 
a physical state corresponding to \a k ) and given by the same expansion coefficients but 
over the eigenstates of the exact Hamiltonian, \a e k ) = Y2e a k{E) \ip e (E)) with a k (E) = 
(ip e {E) \a e k ) . We can therefore understand non- Hermit icity as a consequence of working 
with exact wavefunctions but defined only over part of configuration space relative to the 
underlying exact problem. 

As was the case with H [Eq (|T2l ] that needed to be replaced with H, an operator A 
Hermitian in T-L p h is represented in % by an operator A whose expansion over the biorthogonal 
basis reads 

A = J2mE))A E E'(HE')\. (21) 

EE' 

The relation between A and A is given by j^, [ljj] A = Q~ x l 2 AQ 1 ^ 2 . The time evolution 
operator is a prominent example: U(t) = 'Yl lE e.~ lEt \ip(E)) (tp(E)\ is not unitary in % p h- The 
correct unitary operator in H p h is obtained from U = Q^^UQ 1 ! 2 where U(t) defined by 

U(t) = Y,e- tEt mE))(^E)\ (22) 

E 

is (pseudo) unitary in H. 
D. Weak values 

We are now in a position to formulate the weak value expression for a non-Hermitian 
system. First note that, as described in Sec. 2, we do not need to solve explicitly the full 
)roblem involving the coupling of a non-Hermitian system to a Hermitian one (as eg in Ref. 



16j). Indeed the non-Hermitian system is practically not affected by the weak interaction, 
while the Hermitian one (the meter) simply picks up a phase. This phase - the weak value 
- is therefore the only quantity we need to determine. 
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Let \cti) = J2E a i(E) \ip(E)) be the initial ("preselected") state prior to the weak mea- 
surement of a system observable A and \/3f) = J^E^fi^) IV'(-E')) t ne "postselected" state 
obtained after the subsequent projective measurement. According to the discussion above, 
\cti) and \(3f) represent the physical states in % and the coefficients a; and bf are given by 
formulae analog to Eq ( I20p . The observable A is represented in % by the non-Hermitian 
operator A whose expansion over the biorthogonal basis was given by Eq (l2~Tj) . 

The formula ([3]) expressing the weak value of the observable A becomes 

\ A /w ~ 7^)7 (23) 

or in terms of the standard inner product notation 



I w 



(fif 


A \cti) 


(fit 


Oti) 



(PAS- 1 ]*) me- 1 ]*) ' 



(24) 



Hence when a system that is non-Herrmitian (relative to the standard inner product) inter- 
acts with a weak measurement apparatus measuring the observable A, the meter is shifted 
by a quantity given by Eq (124]) . not by Eq (EJ); Eqs (124"]) and (j3J) obviously coincide when 
the metric is flat (Q is the identity operator). Note that the shift in the weak measurement 
apparatus can in principle be experimentally observed. 

IV. COMPUTATION OF WEAK VALUES IN A MODEL NON-HERMITIAN 
SYSTEM 

A. Model 

We will give examples involving the computation of Eq ( 124|) in a situation well-known in 
atomic physics involving atoms with a single excited electron. In this case the long-range 
field is the familiar Coulomb potential and the reaction zone is about the size of the atomic 
core. The excited electron periodically scatters off the core, the core-electron interaction 
being embodied in the short-range potential. We set up a model with 5 scattering channels: 
the target has a ground state with e\ — and 4 excited states with energies e^i = 2, ...,5. 



The 5x5 scattering matrix K(E) is chosen to have a very strong energy dependence 22] 



in order to have stronger non-diagonal elements of the metric Q (a non-Hermiticity index 
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k can be denned by averaging over the N largest non-diagonal elements of Q where N is 
the dimension of the metric; in this model we have k = 0.061, a non- negligible though 
relatively samll value). The bound state energies E are obtained numerically by enforcing 
the boundary conditions in Eq (fTUj) and then the coefficients Z^E) are retrieved by solving 
the relevant linear system. While the number of bound states is infinite, good numerical 
convergence is obtained by taking about 200 states above and 200 states below the energy 
interval of interest. The metric employed in the numerical computations is thus a 400 x 400 
matrix. 

For the purpose of illustration we will determine the weak value of the radial position of 
the excited electron and the weak value of the energy, assuming in both cases postselection 
can be made to a final state identical to the initial one. We choose an initial state \a(t = 0)), 
that we take to be a Gaussian localized radially very far from the target, at the outer turning 
point of the radial potential for an excited electron (with a mean energy n = 42), with the 
target being in its ground state. Initially \a(t = 0)) is defined on an orthogonal basis of H p h 
but we assume (and verify numerically) that this state can approximately be expanded on 
our chunk of computed eigenstates of H as 

\a(t = 0)) = \F loc (r « r tp )) = £ a(E) \i/)(E)) (25) 

where the a(E) are determined as in Eq ( 1201 . We now proceed to compute weak values. 



B. Weak value of the energy 

We consider a scheme in which a weak measurement of the energy is made at t = tw, 
immediately followed by a projection to a final state. We assume for definiteness it is 
possible to postselect on a state \/3f) identical to the initial state \a(t = 0)), for example 
by considering a weak measurement apparatus consisting in an array of devices placed 
spherically at a radial distance r ps r tp from the atomic core; the weak measurement time tw 
must then correspond to the recurrence time (when the wavepacket relocalizes periodically 
at the turning point ljj]) in the initial scattering channel, here channel 1. 

If non-Hermitian issues are ignored, then the operator H of Eq (|12p would be employed 
for the Hamiltonian, the evolution operator, accounting for the evolution of the system form 
t = to tw would be taken as U(t) = J2e e~ lEt \if}(E)) {%fj{E)\ and the weak value obtained 
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from the usual definition would thus be given by 

(H( , » «_„ (p f \Hm w )) 

mtw))w ~ W\ *M) ( } 

u = „ (a(t = 0)| gj/faQ \a(t = 0)) 

(a(t = 0)\U(t w )\a(t = 0)} ' 1 J 

This quantity is plotted in Fig. 1 (dotted lines) for different possible choices of the measure- 
ment time tw compatible with the system wavepacket radially localized in the neighborhood 
of the measuring apparatus. 

However, since the system is non-Hermitian, Eqs ()26l)-(l27l) should formally be replaced 

by 



(flf 


H\iP(t w )) 


(Pf 


ip(tw)) 



\ HM / W = — 

_ (a(t = 0)\HU(t w )\a(t = 0)) 



(5(t = 0)| U(t w ) \a(t = 0)} 

where H is the non-Hermitian Hamiltonian given by Eq ( !T6l) and U (t) is the corresponding 
evolution operator given by Eq ( 122|) . This quantity is also plotted in Fig 1 (solid line). 

The results shown in Fig. 1 indicate a similar overall behaviour for the two curves, though 
there are substantial differences for several values of the measurement time 23( Therefore 



the replacement of the usual formulae (j25]) - (j2Tjl by Eqs. (T28T) - f|29T) is not purely formal: in 
practical computations the non-Hermitian nature of the system, coupled to a weak measure- 
ment device, must be taken into account in order to compute correctly the expected shift 
in the pointer of the measurement apparatus due to the weak measurement. 

C. Weak value of the momentum 

Another example is the weak value of the momentum postselected to a given position. This 



has become a standard example 18| involving a weak measurement of an observable which 
is incompatible with the postselected one. Here, rather than postselecting to a position \r) 
known with an infinite precision, we employ as above \(3f) = \a(t = 0)) as the postselected 
state, keeping in mind that (r| a(t — 0)) is a wavefunction tightly localized around the 
turning point r tp [Eq ()25l)]. As in the previous example we assume the weak measurement 
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FIG. 1: The weak value of the energy (in atomic units/10~ 5 ) is shown for different measurement 
times compatible with the system wavepacket being in the neighborhood of the weak measurement 
device. The dashed (gray, online red) line represents the usual weak value expression, given by 
Eqs (|26p -(|27p. The solid black line represents the correct expressions Eqs (|28p - (|29p for the weak 
value of the energy in non-Hermitian systems. The time is given in units of the wavepacket period 
(about 1.15 x 10" 11 s). The average energy of the system is 27.89 x 10 -5 au. 

on the preselected state can be made at different times t\y for which the system wavepacket 
relocalizes in the neighborhood of the measuring apparatus. 
The weak value (J24l) becomes 



P{tw)) =A t^ ( 3 °) 



with \ip{tw)) = U{tw) \ a (t — 0)). Note that (P(tw) ) has both a real and a complex 
part: the real part is related to the average velocity field of the system while the complex 





pm w )) 


(Pf 


ii>{t w )) 



part is proportional to the logarithmic derivative of the system wavefunction modulus [19]. 
In principle both the real and the complex parts can be experimentally observed (though 
not jointly). The real part of the weak value (!30|) is plotted in Fig. 2 (black solid line). 
The dashed line is obtained by a straightforward application of Eq fl3]), ie when the non- 
Hermitian character of the system is not taken into account. The two curves nearly overlap, 
which can appear as a little surprising in view of the fact that U(t) is not unitary and 
therefore probability is not conserved. Notwithstanding there are measurement times for 
which the discrepancy between Eq ([3]) and the correct Eq (|24|) is important. 
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FIG. 2: The weak value of the momentum (in atomic units/10 -3 ) is shown for different mea- 
surement times compatible with the system wavepacket being in the neighborhood of the weak 
measurement device. The dashed (gray, online red) line represents the usual weak value expression 
while the solid black line shows the expression valid for non-Hermitian systems. The time is given 
in units of the wavepacket period (about 1.15 x 10 -11 s). 

V. SUMMARY AND CONCLUSION 

We have investigated weak measurements for quantum systems described by a non- 
Hermitian Hamiltonian. The standard definition fl3]) of the weak value - the quantity that 
can in principle be experimentally observed by reading the pointer of a weak meter - does 
not hold in a non-Hermitian framework. The modified expression, given by Eq ( )24"|) was 
derived in this work in the context of bound state scattering theory. The advantage of em- 
ploying this particular instance of non-Hermitian formalism is that its physical meaning is 
devoid of any ambiguity, though the validity of Eq (124"]) holds in general (at least when the 
relevant similarity transform can be denned). 

The results were illustrated numerically in a model system by computing the weak values 
of the energy and of the momentum, with a postselection to a state identical to the initial 
radially localized wavefunction. Overall, the results indicate that even in a non-Hermitian 
in which the non-diagonal elements of the metric are relatively small (the non-Hermiticity 
index was k = 0.061 <C 1), it is important in concrete studies of weak measurements to 
employ the correct (ie non-Hermitian) formulae in order to account appropriately for the 
behaviour of the weak meters. 
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